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Abstract 

We study the physical features of a class of exact solutions for cold 
compact anisotropic stars. The effect of pressure anisotropy on the maxi- 
mum mass and surface redshift is analysed in the Vaidya-Tikekar model. 
It is shown that maximum compactness, redshift and mass increase in the 
presence of anisotropic pressures; numerical values are generated which 
are in agreement with observation. 
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1 Introduction 

The description of very compact astrophysical objects has been a key issue in 
relativistic astrophysics for the past decades. Recent observations suggest that 
there are many compact objects such as X-ray pulsar Her X-1, X-ray burster 
4U 1820-30, miUisecond pulsar SAX J 1808.4-3658, X-ray sources 4U 1728-34, 
PSR 0943-1-10 and RX J185635-3754, whose estimated masses and radii are not 
compatible with the standard neutron star models. The conjecture that quark 
matter might be the true ground state of hadrons [U [5] , inspired many authors 
to describe such stars as strange stars |S1 |3], quark-diquark stars fS], hybrid 
stars [6] and boson/boson- fermion stars [7l [51 [9l flOl [TT] . 
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As densities of such compact objects are normally above nuclear matter 
density, theoretical studies suggest that pressures within such stars are likely 
to be anisotropic, i.e., at the interior of such stars there are two different kinds 
of pressures, viz., the radial pressure and the tangential pressure [12'. Differ- 
ent solutions of Einstein's field equations for anisotropic fluid distribution with 
spheroidal geometry, with varying forms of the energy density, have been ob- 
tained by many workers [131 IHl [HI [S]. So far, the role of pressure anisotropy 
has been extensively studied in the context of high redshift values and stabil- 
ity of compact objects (see for example [T71 [T51 [T^] and references therein). 
Bowers and Liang [5D| have pointed out that anisotropy may also change the 
limiting values of the maximum mass of compact stars. The objective of the 
present work is to investigate the role of pressure anisotropy on the maximum 
masses of compact objects. To this end, in §2, we modify a solution obtained by 
Mukherjee et al [21] to incorporate anisotropy. The class of solutions, capable of 
describing cold compact stars, was obtained by using an ansatz given by Vaidya 
and Tikekar [22j . For physically relevant anisotropic stars, the regularity and 
matching conditions for the solutions are developed. In §3 we discuss the role of 
anisotropy and calculate the maximum possible masses for this class of solutions 
in §4. We conclude by summarizing our results in §5. 



2 Anisotropic model 

We take the line element for a static spherically symmetric cold compact star 
in the standard form 

where 7(r) and ^{r) are the two unknown metric functions. Assuming the 
energy momentum tensor for an anisotropic star in the most general form 

T,j ^ diag (-p, pr, p±, p±) (2) 

the field equations are obtained as 

P = ^ — ;:2 — + — - — ' W 

Pr = ^^-^^f^, (4) 

Ae^^ = 7" + 7'^-7y-^-^-^^, (5) 

where we have set p± — Pr = A. In ([21)-®, p is the energy density, pr is the 
radial pressure, pj_ is the tangential pressure and A is the measure of pressure 
anisotropy in this model. To solve this system we use the ansatz [22^ 

l + VV^ ,2^,„ll 
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Then ^ takes the form 



(1 + A - Xx^)'^xx + Ax*a; + A(A + 1)^' 



2\2 



* = 0. 



(1-X2) 

To solve ([7]) we assume that the form of the anisotropic parameter A is 

aA2(l -a;2) 



(7) 



A = 



i?2(l + A- Aa;2)2^ 



and if wc make a further transformation z = X/{X + ([7]) becomes 

(l-z2)*,, + z*, + (A + 1)* = (8) 

where a = 1 — A/A is a constant. Tliis has the general solution (see ^T] for 
details) 

^ _ r cosp+l)C + <5] _ cosp-l)C + ^] 1 

where, [3 = \/A + 2, = cos^^ z, and A and (5 are constants which can be 
determined from the boundary conditions. The physical parameters in this 
model are then obtained as 
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P = 



Pr 

P± 
A 



i?2(l-z2) 
1 

~i?2(l_^2) 
Pr + A, 



1 



(A + l)(l-z2) 



(A + 1)* 



aX 



iX + l){l-z')-l 



(A + l)2(l-z2)2 



(10) 

(11) 
(12) 
(13) 



which together with ^ and ^ comprise an exact solution to the Einstein field 
equations. Note that 

(1 + A)63 



M{b) = 



2i?2(i + A§; 



(14) 



is the total mass of a star of radius b. 

We impose the following conditions in our model: 

• At the boundary of the star the interior solution should be matched with 
the Schwarzschild exterior solution, i.e.. 



g27(r=6) ^ g-2M(r=fc) 



2M\ 



(15) 



• The radial pressure Pr should vanish at the boundary of the star which 
gives 

III. Cr,^ n 4- X^i 

(16) 



^z{zb) _ (1 + A) 



*(^6) 



2zb 
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where = (A/(A + 1))(1 - b'^/R^). From ® we have 



_ iP^ - 1) 



sin[(/3 - 1)C + 6]- sm[{p + 1)C + S] 



il3 + l) cos[(/3 - 1)C + <5] - (/3 - 1) cos[(/3 + 1)C + S] 

Combining (fT6|) and p7|) we obtain 

r cot^h — tan(/3(^f,) 



tan (5 



1 + root Chtan(/3C6) 



(17) 
(18) 



where 



A(l-2a) + l 



and Cb = cos Zb- 

• Pr > inside the star gives 

1±<J1±R 
^ - 2z ' 

• Using (Unil-dni) we get 

dpr Z(l - z2)2(5f J^r)2 _ (1 _ z^)^'J^) - aXz{l - Z^) 



dp 
dp 



dpr 

dp ^ (1 + A) 



z(l - z2)(l + A) +4z 
(A + l)(l-z2)-2' 
(A + l)(l-z2)+4 



(19) 

(20) 
(21) 



We choose the parameters so that the causahty conditions are not violated, 
^ < 1 in this model. 

The above conditions are imposed for a physically reasonable model. 



3 Physical applications 

It was shown earlier by Sharma et al pS] that the Vaidya-Tikekar model provides 
a simple method of studying systematically the maximum mass problem of 
compact isotropic {a — 0) stars. To see the effect of anisotropy (a 7^ 0) in this 
model, we may adopt the following methods. 

We may choose the isotropic compactness Ui — {M/b)iso and A as input 
parameters and using (fT4|) calculate y — b^/R^. For a given central or surface 
density, (fTO|) can be used to calculate the value of R which then determines the 
radius b = R^ or mass M (from ([U])). The parameter 5 is fixed by choosing a 
specific value of a. Since mass and radius are fixed, this method is not suitable 
to analyze the role of anisotropy on the maximum mass problem. However, 
(PT|) and (PT|) can be utilized to show that for stars with the same masses and 
radii may have different anisotropic compositions if the equations of state are 
modified accordingly. In §3.1, we consider two such examples to show how the 
composition may change in the presence of anisotropy. Note that the variations 
of the slopes of and may correspond to different material compositions 
within the star. These variations are shown in figure [T] and figure [21 
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A = 2, = 7.7 km 


A = 53.34, b = 7.07 km 


a 


Vani 




Mam (Mq) 


a 


Vani 




Mam (Mq) 





0.1450 


0.1686 


0.88 





0.0267 


0.2994 


1.435 


0.2 


0.1281 


0.1530 


0.80 


.05 


0.0269 


0.3002 


1.439 


0.4 


0.0953 


0.1201 


0.63 


0.1 


0.0270 


0.3006 


1.441 


0.5 


0.0690 


0.0909 


0.47 


0.2 


0.0268 


0.2997 


1.437 


0.6 


0.0322 


0.0454 


0.24 


0.4 


0.0239 


0.2858 


1.368 


0.65 


0.0082 


0.0121 


0.06 


0.6 


0.0090 


0.1760 


0.844 



Table 1: Compactness and mass calculated for different anisotropic parameters 
for two different cases discussed in §3.1. 



To see the effect of anisotropy on the compactness, we adopt a different 
approach. Note that (jl8p is a relation between y and S which we will utilize to 
calculate 6 for given values of a and A. For given values of A and Ui, we first 
calculate y using (HH). Substituting these values in (US]), we determine 6 for the 
isotropic case (a = 0). Once 6 is determined, we use this value to calculate yam 
for different a values. We then use the relation 

_ (1 + X)yani 

2(1 + Ay„.) 

to see the effect of anisotropy on the compactness of a star. 
3.1 Numerical results 

Following the method discussed above, we have obtained numerical results show- 
ing the effect of anisotropy on some physically relevant parameters. Two differ- 
ent cases have been studied. 

Case I: We use our earlier data for the pulsar Her X-1 [24] and choose A = 2, 
M — 0.88 Mq, b = 7.7 km so that Ui — 0.1686 and calculate the compactness 
for different anisotropic parameters. 

Case II: We consider the millisecond pulsar SAX J 1808.4-3658 and use the 
results obtained in an earlier work 25j and choose A = 53.34, M ~ 1.435 Mq, 
b = 7.07 km so that = 0.2994. 

We note that the compactness decreases with increasing anisotropy which is 
in agreement with earlier results obtained in [17j . The results are shown in table 
1. The behaviour of the anisotropy factor (A = R^A) in the stellar interior is 
shown in figure ([3]). 

4 Maximum mass and Surface redshift 

In an earlier work |23] , we calculated the maximum mass for a class of isotropic 
stars. Here we follow the same technique to calculate the maximum mass in the 
presence of pressure anisotropy. 
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0.15 

0.1 0.2 0.3 0.4 0.5 0.6 

a 

Figure 1: Variations of ^ at the boundary and at the centre of an anisotropic 
star against a. We took A = 2 and Uj = 0.1686. The soUd line is for (^)r=0) 
the dotted Hne is for (^')r=o, the long dashed Une is for and the 

dashed hne is for {^^)r=b- 




Figure 2: Variations of ^ at the boundary and at the centre of an anisotropic 
star against a. We too^ A = 53.34 and Ui = 0.2994. The soHd line is for 
(^)^=o, the dotted line is for (^)r=o, the long dashed line is for (^)^=(„ 

and the dashed line is for {^^)r=b- 
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. 15 




0.76 0.77 0.78 0.79 0.8 0.81 

z 

Figure 3: Variation of anisotropy factor (A = R^A) against radial parameter z 
for Ui = 0.1686, A = 2 and a = 0.6. 



a 


Umax 


^max 


■^s 1 max 


6 = 10 km 


pb = 2 Pnucl 










M^ax {Mq) 


M„,ax (Mq) 


hmax (km) 





0.0252 


0.3615 


0.9003 


2.45 


2.60 


10.62 


0.2 


0.0285 


0.3738 


0.9910 


2.53 


2.69 


10.62 


0.4 


0.0322 


0.3852 


1.0872 


2.61 


2.77 


10.62 


0.6 


0.0361 


0.3955 


1.1879 


2.68 


2.84 


10.62 


0.7 


0.0383 


0.4003 


1.2398 


2.71 


2.88 


10.61 


0.8 


0.0404 


0.4048 


1.2927 


2.74 


2.91 


10.60 


0.9 


0.0427 


0.4092 


1.3463 


2.77 


2.93 


10.59 


1.0 


0.0450 


0.4132 


1.3998 


2.80 


2.96 


10.58 



Table 2: Maximum compactness {umax), maximum surface redshift {Zs\max) 
and maximum mass {Mmax) for different anisotropic parameters . We have 
considered a star of radius 10 km and surface density equal to twice Pnucl, 
where, pnud = 2.7 x 10^^ gm/cm^. 
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• We assume that < 1 and the value is maximum at the centre. This 
gives 

(1 + A) 



> 



2^/A 



VA+T 



21A+ 1 



4aA2 



A + 1 



(22) 



Combining (|17p and (|22|) , we determine the Umiting value of S for different 
a values for a chosen value of A. 

• Corresponding to the limiting value of 5, ([T5)l can be used to calculate the 
maximum value of y = / R^. 

• From ()14p the compactness of a star in this model is given by 

M{b) (1 + A) 



2(A+i) 



(23) 



Clearly, the maximum value of y corresponds to the maximum compact- 
ness of the configuration. The maximum surface redshift {Zs\max) corre- 
sponding to this value can also be obtained using the following equation 



Z s\max (1 '^'^ani) 



-1/2 



1. 



(24) 



Once the value of maximum compactness is obtained, the maximum mass of 
anisotropic star can be calculated for a given radius or surface density. In [23j 
we observed that for a particular choice (A — 100), the maximum compactness 
for an isotropic star is 0.3615. Keeping the same value of A if we go on increasing 
a we see that the maximum compactness, maximum surface redshift and maxi- 
mum mass all increase with anisotropy. The results are shown in table 2. For a 
close to unity (the maximum value of a in the present model is 1) these values 
are almost 0.4, 1.4 and 2.8Mq, respectively, for a star of radius 10 km. These 
values are similar to the results obtained in [26j . The maximum surface redshift 
obtained by Bondi [27] was 1.352 which is also very close to our values. The 
maximum mass for an isotropic star of radius 10 km was 2.45Mq |23| . which 
increases to 2.8Mq in the presence of anisotropy. 



5 Discussion 

We briefly point out the behaviour of the dynamical variables in this class of 
models. It is clear that the energy density p and the radial pressure Pr are 
decreasing functions from the centre to the boundary of the star. This is also 
true for the anisotropic stellar models of Chaisi and Maharaj [12] and Sharma et 
al [23] who have studied the same spacetime geometry. The tangential pressure 
p± has a more complicated behaviour because it is related to the anisotropy 
factor via p± = pr + A; in addition p± depends on the Gegenbaur function and 
the new variable z rather than the original radial coordinate r. To illustrate 
the behaviour of p± we have generated a plot in Figure 4. It is clear that the 
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Figure 4: Radial pressure (solid curve) and tangential pressure (dashed curve) 
plotted against the radial parameter z for a — 0.4 and A = 2. 

tangential pressure is an increasing function as we approach the centre. This is 
physically acceptable since the conservation of angular momentum during the 
quasi-cquilibrium contraction of a massive body should lead to high values of 
p± in the central regions of the star. 

We have extended a class of solutions describing cold compact stars to incor- 
porate anisotropy. The solutions were then used to see the effect of anisotropy 
on the maximum possible mass and surface redshift parameters of cold compact 
stars. A comparative study of our results with earlier results are given in table 3. 
The anisotropy in the present model vanishes at the centre and reaches the maxi- 
mum value at the surface of the star as shown in the figure 3. Unlike some earlier 
works [18', 19^ , this model has an isotropic counterpart (a ~ 0) which helps to 
compare anisotropic stars with their isotropic counterparts. In this model we 
assumed p_L > pr and have shown that the upper bound on the maximum mass 
increases in the presence of anisotropy. To conclude, our model provides a sim- 
ple method to fix the upper bound on the maximum possible masses for the 
class of compact anisotropic stars described by the Vaidya-Tikekar model. 
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Guven and Murchadha [28j 


0.974 


5.211 


Ivanov 1291 


0.957 


3.842 


Bondi [221 


0.819 


1.352 


Hernandez and N-itnez j26j 


0.800 


1.200 


Present Work 


0.826 


1.400 



Table 3: Maximum compactness {umax) and maximum surface redshift {Zs\max) 
of anisotropic stars in different models. 
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